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Cosmic rays: current status, historical context 
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The ISVHECRI conference series emphasizes the connection between high energy physics and cosmic ray 
physics — the study of elementary particles and nuclei from accelerators in the lab and from space. In this 
introductory paper on cosmic rays, I comment on several current topics in the field while also providing some 
historical context. 
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1. INTRODUCTION 

The first meeting in the ISVHECRI (International 
Symposium on Very High Cosmic Ray Interactions) 
series was in 1980 in Siberia. It is my belief, however, 
that the Bartol meeting in 1978 jlj really set the stage 
for the series by bringing together some key people 
from particle physics, including Carlo Rubbia, David 
Cline, J.D. Bjorken, Gaurang Yodh, Larry Jones and 
Francis Halzen, and from cosmic ray physics, includ- 
ing S. Miyake, C.M.G. Lattes, K. Kamata, Y. Fuji- 
moto, Michael Hillas, K. Niu and W. Vernon Jones. 
Many of the same people were together at the second 
conference in the series in 1982 which started in La 
Paz and ended in Rio de Janeiro [2j. Gaurang Yodh 
includes a color print of the photo from the Bartol 
conference in his Hess Lecture presented at Pune in 
2005 jsj , along with several other historic photos that 
illustrate his comprehensive review of cosmic rays. 

My charge at this conference is to provide introduc- 
tory comments on the status of cosmic rays. The orga- 
nization of the paper is conventional, starting with in- 
clusive measurements of atmospheric muons and neu- 
trinos, followed by direct measurements of primary 
cosmic rays at the top of the atmosphere and finally 
extensive air shower measurements in the PeV region 
and above where the intensity is too low for direct ob- 
servations. Most of these topics are covered in detail 
in papers presented during the conference. I therefore 
focus on a few points of interest without attempting 
a comprehensive review. 



2. MUONS AND NEUTRINOS 

An approximate analytic expression for the inten- 
sity of atmospheric leptons is 



ME L ) = <t> N (E t ) x 
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where 4>n(E l ) — dN/d\n(E L ) is the primary spec- 
trum of nucleons (N) evaluated at the energy of the 



lepton (t). The forms are the same for i — v and 
l = /-i with different numerical values for the constants 
Aj L etc. [1, H[. The three terms correspond respec- 
tively to contributions from decay of pions, kaons and 
charm Q. The approximations assume that the pri- 
mary spectrum of nucleons can be well approximated 
by a power law, that hadronic cross sections are en- 
ergy independent, and that inclusive cross sections for 
hadro-production obey Feynman scaling. The devel- 
opment of this analytical approach dates back 50 years 
to the early papers of Zatsepin 7] and others. Al- 
though detailed simulations are needed in the end, the 
analytic approximations provide quantitative insight 
into the physics of atmospheric leptons. 

The classic paper of Barrett et al. [8] deals with 
measurements of muons in deep detectors that use 
the overburden to select muons which had high energy 
at production in the atmosphere. They describe how 
the angular dependence of the muon intensity deep 
underground and its dependence on temperature in 
the stratosphere depend on the relative contributions 
of 7T^~ and K ± to the muons. Both effects arise from 
the competition between decay and re-interaction of 
the parent hadrons in the upper atmosphere. The 
expression for the critical energy for meson decay is 



E, 



critical 
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(2) 



where 0* is the local zenith angle at lepton production 
taking account of the curvature of the Earth |5( . The 
parameters are e w « 115 GeV for charged pions as 
compared to ejf as 850 for charged kaons. 

2.1. Temperature Dependence 

In an isothermal approximation, the density of the 
atmosphere is described by an exponential with a scale 
height of Iiq as 6.4 km, where the numerical value is 
applicable to the stratosphere where most high-energy 
muons and neutrinos originate. This numerical value 
is used to compute the values of Cj in Eq. [2] that ap- 
pear in the denominator of Eq. [TJ From the ideal 
gas equation relating density and pressure, one finds 
ho = RT. Since e oc ho, this relation gives rise to 
a seasonal variation in the rate of muons in a deep 
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underground detector through the dependence of e in 
Eq. Q]on ho- Rates can also reflect sudden changes in 
the upper atmosphere [9(. 

The magnitude of the correlation with temperature 
depends on the muon energy at production and hence 
on the depth of the detector. Because ex > Ctt, 
the kaon contribution has a weaker correlation with 
temperature than the pion contribution, which means 
that the magnitude of the correlation with temper- 
ature depends on the kaon to pion ratio. This de- 
pendence has been discussed recently by the MINOS 
Collaboration in connection with their measurements 
of cosmic-ray muons fioj in the far detector at Soudan, 
where the minimum muon energy to reach the detec- 
tor is approximately 1 TeV. In IceCube the seasonal 
variation in the muon rate is approximately ±8% [Tlj . 
This is an important feature of the background that 
must be accounted for in interpreting the data. 

Rates of atmospheric neutrinos are also expected to 
vary with temperature [T3 |. IceCube is large enough 
so that the rates of TeV neutrinos may be sufficiently 
high to actually measure their seasonal variation Q. 
Since the neutrinos are detected primarily as upward 
moving ^-induced muons (using the earth as a filter 
against the downward muon background), the magni- 
tude and phase of their variation will depend on the 
part of the sky from which they originate. 

Because the lifetime of charmed hadrons is so short, 
their decay products are "prompt." The denomina- 
tor of the charm contribution in Eq. [T] is unity for 
E L < 10 PeV. As a consequence, the prompt muons 
and neutrinos are isotropic and have a harder spec- 
trum than decay products of pions and kaons in the 
energy range above a TeV. These two features have 
long been used to search for aprompt contribution to 
atmospheric muons. In Ref. [6] it is pointed out that 
the lack of seasonal variation is also a signature of the 
charm channel. 

2.2. Muon Charge Ratio 

Both MINOS 13] and OPERA QJ] are sufficiently 
deep so that they measure the charge ratio of muons 
that had energies in the TeV range at production in 
the atmosphere. The measured increase in the charge 
ratio from 1.27 below 100 GeV to 1.37 above a TeV 
reflects the enhanced importance of kaons at high en- 
ergy. In particular, associated production, 

p + air -)• A + K + + X, (3) 

is an important factor in making the +/— charge ratio 
larger for kaons than for pions. A measurement of 
the charge ratio by CMS at somewhat lower energy 
was reported at this conference [HI . The results from 
MINOS on temperature dependence as well as charge 
ratio were summarized here by Schreiner [llij . 



The strong asymmetry of K + compared to K~ in 
the forward fragmentation region (which dominates 
production of secondaries from a steep spectrum) re- 
flects the fact that the A has constituents in common 
with the incident proton, while there is no correspond- 
ing channel for production of K~ . The charmed ana- 
log of Eq. [3] is also highly asymmetric pjj • The pos- 
sible consequence for prompt atmospheric leptons is a 
subject of current interest. 

Analytic expressions for charge separated intensities 
(/i + , and v^) can be constructed in parallel 
with Eq. [I] as was worked out for pions in Ref. [l8| ■ 
The charge-sign dependence starts with the separa- 
tion of neutrons from protons in the primary cosmic- 
ray beam and tracks the charge dependences of meson 
production and decay 

The new data on temperature dependence and 
charge sign dependence of atmospheric muons call for 
a global fit to all the relevant parameters. These in- 
clude all the spectrum- weighted moments, Zj Ll the 
attenuation lengths and also the parameters that re- 
flect the initial charge-sign dependence imparted by 
the fraction of primary nuclei with Z > 1 and how 
this changes with energy. 



3. GALACTIC COSMIC RAYS 

3.1. Direct measurements of the primary 
spectrum 

Measurements by ATIC pj| and CREAM Hj pre- 
sented by Seo at this meeting (2l[ are providing a 
new and more complex picture of primary composi- 
tion in the multi-TeV energy region than was gen- 
erally assumed in the past. In particular, spectra 
seem to harden somewhat above 200 GeV/nucleon, 
and there is a hint that the spectrum of helium is 
somewhat harder than that of protons between 10 and 
100 TeV HJ. 

Figure Q] shows the situation a few years ago. Up to 
100 GeV there are rather precise spectrometer mea- 
surements of protons and helium from BESS (22j, 
AMS [H, CAPRICE [H and others. The measure- 
ments of BESS and AMS agreed well with each other 
in normalization as well as spectrum, with proton in- 
tensities somewhat higher than CAPRICE and earlier 
measurements. The line in Fig. [T] is an extrapolation 
to high energy of the sum of the individual spectra 
assuming that all are described by a power law with 
integral spectral index 7 = 1.7. Normalization of pro- 
tons and helium are to BESS and AMS. Normaliza- 
tion of other nuclei are obtained from HEAO [25| and 
other measurements as summarized in (26j. The ex- 
trapolation comes in on the low side of the indirect 
measurements with air shower experiments around a 
PeV, just below the knee. 
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Figure 1: Summary of measurements before 2005. 



The plot also shows the early measurements of 
JACEE and RUNJOB in the region from 10 to several 
hundred TeV. There was a significant discrepancy be- 
tween JACEE and RUNJOB for the flux of helium. 
ATIC-2 and CREAM seem to confirm the higher He 
intensity of JACEE. Adding together the CREAM 
measurements of protons, helium and heavier nuclei 
gives a higher all-particle intensity around 100 TeV 
in better agreement with the TIBET air shower mea- 
surement. The harder spectrum reported by CREAM 
in the range 10 - 100 TeV/nucleus when extrapolated 
to the knee region leads to this improved agreement. 

3.2. Knee of the spectrum 

Fifty years ago Peters [27( suggested that the steep- 
ening in the knee region may be associated with the 
upper end of the spectrum of particles from the dom- 
inant source of cosmic rays, which could be powered 
by supernovas. He explained the relatively smooth be- 
havior of the cutoff as a consequence of different pri- 
mary nuclei having cutoffs at different energy per par- 
ticle under the assumption that the underlying effect 
depends on magnetic rigidity. The definition of rigid- 
ity, R(Z) = Pc/Ze, provides the scaling relation be- 
tween total energy (E m Pc) and the charge Ze of nu- 
clei moving in a given configuration of magnetic fields. 
If the "cutoff" of the spectrum occurs at a character- 
istic rigidity, R c , then protons bend first at E c = eR c 
followed by helium at 2ei? c , carbon at 6eR c up to iron 
at 26ei? c . This sequence is sometimes referred to as a 
"Peters cycle" . The KASCADE Experiment [H pro- 
vided the first observation of this sequence of changes 
in the knee region, with protons steepening first then 
helium and then heavier nuclei. The relative abun- 
dances of various nuclei inferred from the KASCADE 
data depends on the interaction model used for sim- 
ulation necessary to unfold the spectra. However, he- 



lium appears more abundant than protons above 1 
PeV whether SIBYLL or QGSjet is used for the un- 
folding. This feature is qualitatively consistent with 
an extrapolation of the direct measurements assuming 
a harder spectrum for helium than for protons. 
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Figure 2: Two component concept of Peters (1960). 

Already at the time of Peters' paper the spectrum 
was known to continue beyond the knee. He suggested 
a second population with a significantly smaller total 
number of particles but with a slightly harder spec- 
trum that continues to higher energy, as shown in 
Fig. E3 from his paper. We now know, however, that 
the spectrum continues for three orders of magnitude 
after the knee with a steeper spectrum than below the 
knee. Not until the ankle does the spectrum harden 
in the way expected by Peters. At least one addi- 
tional Peters cycle with a cutoff at higher energy is 
needed to fill in the gap between the knee and the 
ankle. Hillas [II] calls this "Component B" . 

3.3. Component B 

The question of the cosmic-ray spectrum above the 
knee was extensively discussed at the two meetings on 
"Physics at the end of the Galactic cosmic ray spec- 
trum" in Aspen in 2005 and 2007. Figure ED illus- 
trates the discussion with schematic outlines of three 
separate components superimposed on the plot of the 
spectrum from the review of Nagano and Watson [3(| ■ 
In this scheme, the second knee is a cutoff associated 
with a second population of Galactic sources of un- 
known origin, and the extra-galactic component dom- 
inates only above the ankle. It is interesting that the 
new spectrum of KASCADE-Grande presented at this 
conference shows a concavity just above 10 PeV 
that could be interpreted as the onset of Component 
B. 

Another point is to estimate the power required to 
produce the B-component. The desired result is ob- 
tained by integrating the spectrum of the B compo- 



C113 



4 



XVI International Symposium on Very High Energy Cosmic Ray Interactions 
ISVHECRI 2010, Batavia, IL, USA (28 June 2 July 2010) 



M. Nagano and A, A, Watson: Ultra hiotvenergy cosmic rays 



E 

> 



□ Akeno lkm ! 

* Tibet 
A Rynjob 

* Proton Satellite 

* JACEE 




O AGASA 

A Haveron" 

ft Yokutsk 

* Stereo Flys Eye 



log(ENERG¥ In eV) 

Figure 3: Schematic picture showing three populations of 
cosmic rays: A is the dominant galactic component, 
perhaps accelerated directly by supernova remnants; B is 
a higher-energy galactic population of uncertain origin; 
£ — Q is the extra-galactic component that dominates 
above the ankle around 3 EcV. 



nent weighted by an estimate of the lifetime of high- 
energy cosmic rays in the galaxy. The expression 

gives the total power in the sources to supply this 
component. Taking 



2 x 10 7 yrs£; 



-0.33 



v Ga , 



laxy 



= 3 x 10 66 cm 



(5) 
(6) 



and a differential source spectral index below the cut- 
off of 7 + 1 = —2.1, leads to a power requirement of 
~ 2 x 10 39 erg/s. 

3.4. Acceleration mechanisms 



a larger velocity difference than the canonical factor 
of 4 at the sub-shock. In diffusive shock acceleration, 
the integral spectral index 7 is related to the ratio 
of velocity of the upstream flow into the acceleration 
region (u\) to the downstream velocity (112) by 



7 = 



u 1 /u 2 



1 



(7) 



In the non- linear case, Ux(d) is a function of distance 
upstream of the shock, with u\(d)/u2 increasing as d 
increases. Higher energy particles can diffuse further 
upstream on average before reversin g di rection; hence 
the concave nature of the spectrum [34| . 

Non-linear diffusive shock acceleration is contrasted 
with the test-particle approximation in which the ef- 
fect of the accelerated particles on the plasma flowing 
through the shock is neglected. Then, for a strong 
shock, ui/u% = 4 and 7 = 1, which produces a spec- 
trum with equal energy content per logarithmic in- 
terval on energy. In the non-linear case, most of the 
energy of the accelerated spectrum is near E max . An 
interesting idea that uses this feature is the proposal 
by Ptuskin &: Zirakashvili 35]. At each epoch in the 
evolution of a supernova remnant, most of the power 
goes into particles at the highest possible energy. As 
the remnant evolves, E max decreases and the highest 
energy particles can then escape upstream (that is, to 
the outside). The spectrum produced by a supernova 
over its lifetime is then the sum of the contribution 
from each epoch. 

A full picture of the sources, acceleration and propa- 
gation of galactic cosmic rays does not yet exist. What 
seems likely to me, however, is that in reality there 
are many sources with different spectra and capable 
of accelerating particles to different maximum ener- 
gies. The complex structure of the spectrum indi- 
cated by the CREAM results below 100 TeV and by 
the KASCADE-Grande results in the air shower range 
may reflect such underlying complexity. 



Another topic featured at the Aspen meetings was 
a discussion of acceleration mechanisms. The classic 
estimate of the upper limiting energy of acceleration 
by shocks driven into the interstellar medium by ex- 
panding supernova remnants is E mlix ~ 100 TeV for 
protons [32J. Two aspects of non-linear shock accel- 
eration are important for building a synoptic picture 
of Galactic cosmic rays. One is the amplification of 
magnetic fields near the shock [33|], which increases 
-Emax oc ZeBR. This allows diffusive acceleration 
driven by supernova remnants to accelerate particles 
to a maximum rigidity one PV or higher, as would be 
expected if the knee is related to the upper limit of 
the dominant population of Galactic cosmic rays. 

The other important effect is the concavity intro- 
duced by the cosmic-ray pressure of the highest en- 
ergy particles on the plasma upstream of the shock. 
The result is that higher energy particles experience 



4. EXTRAGALACTIC COSMIC RAYS 

A well-known alternative to the picture described 
in the previous section and summarized in Fig. [3] is 
that of Berezinsky and collaborators [3g|. They ex- 
plain the ankle as a consequence of energy loss by 
protons due to e ± pair production on the cosmic back- 
ground radiation. Such a model requires that the 
extra-galactic component consists primarily of pro- 
tons and that there is a cosmological distribution of 
sources capable of accelerating particles to energies 
above 100 EeV. In this picture the extra-galactic pop- 
ulation dominates down to ~ 300 PeV and a galactic 
B component is not needed. A signature would be a 
transition in the composition from heavy nuclei at the 
end of the galactic population to protons that domi- 
nate the extra-galactic component in this model. 
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4.1. Composition 

Results from HiRes [37[ and HiRes-MIA [H do in- 
dicate a transition from heavy nuclei below 300 PeV 
to protons in the EeV range and above. But the ques- 
tion is unsettled because Auger results 39] indicate 
on the contrary that the composition above 10 EeV is 
mostly heavy nuclei. 

The two contrasting pictures of the transition from 
galactic to extra-galactic cosmic rays and their rela- 
tion to composition have been discussed in a series of 
papers by Allard et al. [4(J. In the case of a mixed 
composition at the source, even if there is a preponder- 
ance of heavy nuclei, a significant fraction of protons 
is expected at Earth as a result of photo-nuclear disin- 
tegration during propagation through the microwave 
background, provided that the sources have a cosmo- 
logical distribution and accelerate particles to suffi- 
ciently high energy. In any case, this picture would 
also be difficult to reconcile with the Auger results 
on composition because there would be large fluctua- 
tions in the Xmax distribution as a consequence of the 
mixed composition after propagation to Earth (pro- 
tons and iron). 

There is considerable experimental activity to ex- 
plore the transition region. Low-energy extensions of 
Auger are now in operation (4l| . and the low-energy 
extension of the telescope array is under construc- 
tion E|L Meanwhile, KASCADE-Grande jp, Ice- 
Cube [43| and TUNKA Q are reaching the transi- 
tion region from below. Concerning the discrepancy 
in composition between HiRes and Auger, it may be 
relevant to note that the selection of events is some- 
what different in the two cases. HiRes reconstructs 
event trajectories using two Fly's Eyes in stereo, while 
the Auger reconstruction uses a hybrid technique that 
involves times of hits in one or more surface detec- 
tors along with the fluorescence signal. Some insight 
on the differences may come from the Telescope Ar- 
ray (45[, which is also a hybrid detector. 



4.2. Two generic models 

One possibility is a machine in which the magnetic 
fields required for acceleration are strong enough and 
in an appropriate configuration so that accelerated 
protons and nuclei lose energy in the ambient radi- 
ation fields by photo production and photo disinte- 
gration respectively before they escape from the sys- 
tem. The acceleration time scale has to be short or 
comparable to the photo-interaction rate. Then if the 
escape time for neutrons is short compared to their 
photo-interaction time scale they will escape and con- 
tribute to the population of cosmic rays. In this class 
of models the extra-galactic cosmic rays would be pro- 
tons from decay of the escaping neutrons. 

In addition, neutrinos produced in the same sources 
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Figure 4: Summary of primary spectrum measurements 
from Linsley's review [511 ]. 



in reaction chains like 



P + 7 



(8) 



will then have a spectrum and intensity related by 
kinematics to that of the cosmic-rays. This class of 
models would be a realization of the Waxman-Bahcall 
limit [46]. Ahlcrs et al. [13] point out that the limits 
from IceCube already constrain models of this type. 

A contrasting possibility is that the extra-galactic 
cosmic rays could be powered by shocks driven into 
an external environment. In this case, whatever ma- 
terial is available would be accelerated, in analogy 
with galactic cosmic rays accelerated in shocks driven 
by supernova explosions. External shocks driven by 
gamma-ray bursts or jets of active galaxies [48| could 
be realizations of this class of models. The composi- 
tion could be mixed and there would be no particu- 
lar reason to expect a high intensity of neutrinos or 
gamma-rays from interactions of the accelerated par- 
ticles. 



4.3. The end of the cosmic-ray spectrum 

The experimental discovery of extra-galactic cosmic 
rays can be attributed to the Volcano Ranch experi- 
ment, built by John Linsley and Livio Scarsi starting 
in 1958 [49| . The array was spread over an area of 
approximately 8 km 2 . In 1962 they detected an event 
estimated to have an energy of 10 20 eV (with large un- 
certainty), too high to be contained in the galaxy joOj . 
In his 1978 review of ultra-high energy cosmic rays, 
Linsley [5l[ combined data from Volcano Ranch and 
Haverah Park 52] in the U.K. into a plot that shows 
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Figure 5: HiRes and Auger spectrum measurements from 
Ref. H. 



the ankle of the spectrum, which he interpreted as the 
onset of extragalactic cosmic rays. Linsley's figure is 
reproduced here as Figure |4j The data clearly show 
the ankle but no sign of the suppression expected from 
photo-production interactions in the microwave back- 
ground if the highest energy particles come from a 
cosmological distribution of sources. 

Linsley's review came at an important time in the 
history of measurements of ultra-high-energy cosmic 
rays. The Fly's Eye detector was under construction 
at Dugway Proving Ground in Utah after the tech- 
nique had been demonstrated for the first time by us- 
ing three new Fly's Eye receivers at the Volcano Ranch 
array to see the fluorescence light from air showers in 
coincidence with Linsley's ground array. (53j . 

It was another thirty years until HiRes, using 
the Fly's Ey e technique, showed a drop in the 
spectru m l54l| consistent with that predicted by 
Greisen [55j and Zatsepin & Kuz'min [56J. The HiRes 
result is shown in Fig. [5] along jvith the latest spec- 
trum from the Auger detector 



571 ] . The two spectra 



agree well with each other within the uncertainties of 
the energy calibration, which are of the order of 20%. 
Moving one of the sets of points along the direction of 
the slanted arrow brings the two measurements into 
agreement. In fact this point illustrates the possibility 
noted already by Linsley of using the GZK feature to 
calibrate the energy assignment. 



4.4. Anisotropy? 

Linsley also pointed out that protons with ener- 
gies of the order of 100 EeV would have sufficiently 
high rigidity to be likely to point back to their sources 
within ~ 100 Mpc. In this connection it is interest- 
ing to compare the sky plot in Linsley's review with 
that of the Auger experiment. Linsley shows the di- 
rections of 16 events with e nerg y greater than 50 EeV 
measured by Haverah Park [58| ■ He points out that a 



significant fraction of the events come in from direc- 
tions at large angle to the galactic plane, which argues 
in favor of an extragalactic origin if the particles are 
protons with high rigidity. The first Auger sky plot 
published in 2007 [59| showed a rather strong correla- 
tion of events with E > 57 EeV with the large scale 
distribution of matter in the nearby universe as traced 
by AGNs. Since then the significance of the correla- 
tion has decreased [6(j, but it is still suggestive of a 
correlation of the kind that might be expected if the 
sources are correlated with the large-scale distribution 
of matter [61j | and if the particles in this energy range 
are mostly protons. 

In this connection, it is important to note that 
the measurement showing heavy composition extends 
only up to 40 EeV and does not yet overlap with the 
event sample used to look for anisotropy. In a re- 
cent paper, Calvez, Kusenko and Nagataki |62j sug- 
gest that GRBs in the Milky Way on a 10 5 year time 
scale could contribute to the pool of cosmic rays above 
10 EeV. This would be a kind of second "B compo- 
nent" that would explain the heavy composition of 
Auger as the end of the highest energy galactic Peters 
cycle. The onset of the true extragalactic component 
would then be around 50 EeV. The challenge is to find 
ways to increase the statistics and precision of mea- 
surements of cosmic rays up to 100 EeV despite the 
extremely low intensity of such particles. 

4.5. End of the spectrum revisited 

Finally, it is salutary to remember the Hillas 
plot [63] . In his iconic diagram, Hillas places potential 
sources in a space defined by their sizes and magnetic 
fields. A remarkable feature of the plot is that the few 
sources capable of accelerating particles to 100 EeV 
are not far above the line. A realistic possibility is, 
therefore, that the steepening of the spectrum may to 
some extent reflect a cutoff at the accelerators as well 
as (or instead of) the effect of energy losses to inter- 
actions with photons of the cosmic microwave back- 
ground. 
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